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Abstract 



The dielectric dispersion, dielectrophoretic (DEP) and electrorotational 
(ER) spectra of spheroidal biological cell suspensions with an intrinsic dis- 
persion in the constituent dielectric constants are investigated. By means of 
the spectral representation method, we express analytically the characteristic 
frequencies and dispersion strengths both for the effective dielectric constant 
and the Clausius-Mossotti factor (CMF). We identify four and six charac- 
teristic frequencies for the effective dielectric spectra and CMF respectively, 
all of them being dependent on the depolarization factor (or the cell shape). 
The analytical results allow us to examine the effects of the cell shape, the 
dispersion strength and the intrinsic frequency on the dielectric dispersion, 
DEP and ER spectra. Furthermore, we include the local-field effects due to 
the mutual interactions between cells in a dense suspension, and study the 
dependence of co-field or anti-field dispersion peaks on the volume fractions. 
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I. INTRODUCTION 



The polarization of biological cells has a wide range of practical applications like manip- 
ulation, trapping or separation of biological cells [1,2]. Dielectric spectroscopy [3], dielec- 
trophoresis (DEP) [4] and electrorotation (ER) [5] offer a unique capabihty of monitoring 
the dielectric properties of dispersions of colloids and biological cells. Under the action 
of external fields, these particles exhibit rich fluid-dynamic behaviors as well as various 
dielectric responses. Hence, it is of importance to investigate their frequency-dependent 
responses to ac electric fields, which yields valuable information on the structural (Maxwell- 
Wagner) polarization effects. The polarization is characterized by a variety of characteristic 
frequency-dependent changes known as the dielectric dispersion, whose spectra arc helpful 
to analyze the inhomogeneous systems, including biological cell suspensions and tissues [6]. 
Many factors exert influences on the effective dielectric behavior of the system such as ori- 
entation of dipoles, surface conductance and the cell shape [7]. However, some factors can 
be dominant at certain ranges of frequencies (for example, the experimental data revealed 
that the low frequency subdispersions were dependent on the cell shape [8]). 

The DEP is used to describe the motion of the particles caused by the dielectric polar- 
ization effects in non- uniform electric field [9]. The DEP force drives the particles towards 
high intensity (positive DEP) or towards regions of minimal field intensity (negative DEP). 
Because of the interaction (i.e., DEP force) between the induced dipole and the external 
electric field, the particles can be levitated in the medium [10]. 

The ER behavior is caused by the existence of a phase difference between the field-induced 
dipole moment and the external rotating field, which results in a torque to cause the particle 
to rotate. In the dilute limit, the ER of individual cell can be predicted by ignoring the 
mutual interaction between the cells, and hence may be considered as an isolated particle 
in rotation [11]. 

In this paper, we will apply the spectral representation theory [12] to investigate the ac 
electrokinetic behavior including the dielectric dispersion, DEP and ER spectra of biological 
cell suspensions. The object of the present investigation is threefold. First, it is instructive to 
consider the effect of the intrinsic dielectric dispersion. Actually, such an intrinsic dielectric 
dispersion often occurs due to the surface conductivity [13,14] or the inhomogeneous struc- 
ture such as the coated shells of biological cells or colloidal suspensions [15]. Secondly, many 
cells exist in the form of the nonspherical shape such as the fission yeast cell [8] and the red 
blood cell [16], and hence the effect of cell shape needs to be considered [17,18]. Thirdly, the 
volume fractions of the biological cell suspensions are generally not in the dilute limit and 
can even exceed 0.1 [6,8,19], thus the local field effect due to the average electrostatic inter- 
action between cells must be taken into account [20]. In view of this, the effective dielectric 
dispersion spectra will be studied based on Maxwell-Garnett type approximation [21], which 
involves an exact calculation of the field induced in the host medium by a single ellipsoidal 
cell and an approximate treatment of its distortion by the electrostatic interaction between 
different cell suspensions. Moreover, the Clausius-Mossotti factor (CMF), which determines 
the DEP and ER spectra, will be modified by replacing the dielectric constant of the host 
medium with the effective one. 

Being beyond this work, we noticed an alternative paper [22], in which the authors con- 
sidered the electrorotation and levitation of spherical cells or colloidal particles with/without 
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the membrane-covered shells in the dilute limit. However, our investigation can be valid not 
only for nonsphcrical cell suspensions without mcmbrane-covcrcd shells, but also for non- 
dilute limit suspensions. Furthermore, we adopt the spectral representation method [12], 
which offers the advantage of the separation of material parameters from the geometric 
information, to simplify the derivation of the analytical expressions for the characteristic 
frequencies and dispersion strengths of the effective dielectric permittivity and the CMF, 
respectively. 



II. FORMALISM 

We consider a composite system in which biological cells of dielectric constant ei with 
the volume fraction p are dispersed in an isotropic host medium with dilectric constant 
e'i — €2 + 02I (jStt/) with j = V— 1, / being the frequency of the applied field. e\ is assumed 
to exhibit an intrinsic dielectric dispersion, i.e., 

= + TTTTT^ ^ i^' 

where ei (cxi) is the limiting high-frequency (low- frequency) dielectric constant (conductiv- 
ity), Aei represents the dielectric dispersion strength with a characteristic frequency j\. 
Biological cells arc further assumed to be spheroidal in shape. Such an assumption is sup- 
ported by Bohrcn et al [23] who suggested that dielectric dispersion spectra of particles of 
arbitrary shape can be approximated by these of spheroidal particles. 

Then we will adopt the spectral representation to investigate the ac electrokinetic be- 
havior such as the dielectric dispersion, dielectrophoretic and electrorotational spectra of an 
inhomogeneous system in which spheroidal biological cells are randomly distributed. 



A. Dielectric dispersion spectrum 



Generally, for cell suspensions of arbitrary shape, the spectral representation can only 
be solved numerically. However, for randomly-oriented spheroidal cell suspensions, within 
the mean- field theory, the effective complex dielectric constant 4 can be written as [21] 

~ ^ 1 + P [6.(1 - L,) + 2Ky{\ - L^y)\ 

where = (ei — e2)/{3 [e2 + L^^ei — €2)]} {k = z,xy) is called the Clausius-Mossotti factor 
(CMF), and [L^y = (1 — Lz)/2] are the depolarization factors along the z {x or y) axis of 
the spheroids. These depolarization factors depend on the aspect ratio q = c/a [21], where a 
(=6), c are the semi-axes of a spheroid along the Cartesian axes. For the prolate-spheroidal 
cell {q > 1), we have < < 1/3; while for the oblate-spheroidal one {q < 1), we have 
1/3 < Lz < 1. Actually, once a g is given, can be obtained uniquely and thus used to 
indicate the shape of spheroidal particles. 

By invoking the spectral representation and introducing the dimensionless parameter 
s = e2/(e2 — Ci), we rewrite Eq.(2) as 
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Ce = £2 ( 1 



Wi W2 



S — Xi S — X2 



where the poles xi and X2 are given as 
1 
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3-2p + 3L,±J{3-2p + 3L,f - 72(1 - p)L,(l - L, 



Correspondingly, the residues Wi and W2 have 



p(l + SL^ - 6x1) p(l + - 60:2) 

— ~ ^ ) *^2 — 



Q{X2 - Xi) 



6(xi - X2) 



(3) 



(4) 



(5) 



It is easy to check that Wi and W2 satisfy the sum rule Wi + W2 = P- 

In order to obtain the analytic expressions for the dielectric dispersion strengths and the 
characteristic frequencies of the dielectric dispersion spectra, we then introduce two contrast 
parameters s = 62/(62 — ei) and t — a2/{cr2 — cri), firstly defined by Lei et al [7], and derive 
an equality for biological cells possessing an intrinsic dispersion. 



where 
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with 



B{x) = 



s-t 



fc2 — fl 



{s-x){t-x) \ /i/c2 J {s-x)x /a 



/1/2 



fc2{x) - /cl(x)' 
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After simple manipulations, we can express Eq.(3) as. 



with 



Wi W2 \ / Wi W2 
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Aeel = W,A{x,) [ 7^ _ e2 ) , /el = U{x,), (13) 

Aee2 = (^^^^^^ - ^2) , fe2 = Uix,). (14) 

Similarly, Ae^i and /ei for i = 3, 4 can be easily got by our replacing Wi and Xi with VF2 
and X2 from above two equations. 

Thus, within the spectral representation, for spheroidal suspensions with one intrinsic 
dispersion, we successfully obtain four characteristic freqencies f^i (i = 1,4) with the dis- 
persion strengths Ae^i in terms of the geometric parameters (L^ and p) as well as physical 
parameters (s, t, Aei, (J2 and €2). 

In Fig.l, we numerically calculate four characteristic frequencies and dispersion strengths 
against Lz- As is evident from the figure, the characteristic freqencies and the corresponding 
dispersion strengths are strongly dependent on the shape of cell suspensions and exhibit 
nonmonotonical behavior with increasing L^. For these four characteristic frequencies, two 
of them are located at frequencies higher than lO^Hz, while another two are at frequencies 
lower than lO^Hz. Such a separation results from the fact that we take into account one 
intrinsic dispersion in the dielectric cell suspensions. Furthermore, due to the nonspherical 
shape of cell suspensions, there are two characteristic frequencies in the high-frequency 
(or low- frequency) region. When the shape is far from spherical (say — > or 1), the 
difference between these two frequencies becomes large. Conversely, for spherical inclusions 
{Lz = 1/3), only two dispersion strengths are nonzero, leading to two-main steps in the 
dielectric dispersion spectra, as will be shown in the next section. 



B. dielectrophoretic and electrorotational spectra 

The spectral form of Eq. (6) also allows us to investigate the dielectrophoretic and elec- 
trorotational behavior. 

It is known that the time average DEP force F{f) and ER torque r(/) cxccrted on 
the spheroidal cell suspensions with ^-orientation parallel to the electric field are expressed 
as [24,25] 

F{f) = 1.5v~e2Re[bz{ei, £2)]! V ^lms)l r(/) = -3v~e2 Im[h{ei, ~e2)]El, (15) 

where v is the volume of spheroidal cell, 'E(^rms) is the root-mean square magnitude of the 
imposed ac electric field and Re[- ■ ■], Im[- ■ ■] represent the real and imaginary parts of bz, 
which can be written under the spectral representation as 

bz{ei, €2) = J • , ./^ , . = ~ ■ - (16) 
3 €2 + L^{€i -€2) 3 s- Lz 

Eq.(16) is independent of the volume fraction p and thus is valid only in the dilute limit [22]. 

However, for non-dihite volume fractions, we must consider the local field effect due to 
mutual interaction between the spheroidal cells and modify Eq.(16) as. 
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Substituting Eq.(2) into (17), we have 



;i-p)L,(l-L,) 



653 _ (3 + 4p + gi^ _ QpL,)P + {p + 6L, - 3pL,)s + 3(p - 1)L2(1 - L,) 
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+ 
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where Si {i — 1,2, 3) are three roots of polynomial equations 



3 3 + 4p + 9L^ 
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6pL^ 2 , P + QL^ 
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while the residues Vi, V2 and V3 can be obtained from 

Vi + V2 + Vs^ 

Vi{s2 + S3) + ^2(53 + Si) + Vsisi + S2) = 



3pL (p-l)L^(l-L,) 



0, 



(18) 



(19) 



(p-l)(l + L,) 



(20) 



V1S2S3 + V2S1S3 + V3S1S2 = -(p - l)Lz(l - L^). 



Note that these residues again satisfy the sum rule Vi + V2 + V3 = — (1 
Introducing Eq.(6) into Eq.(18) leads to 
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(21) 



The normahzed DEP force and normahzed ER torque can be defined as 



NDEP^Re [6.(61,62)] = 

i=l 



V, 



+ 
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and 
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j^\VA{s,)f/U{s,) , VB{s,)f/f,2{s, 



1=1 



l + [///cl(s, 



1 + [f/fc2{s^)r' 



(22) 



(23) 



Based upon the spectral representation, the expressions for the dispersion strengths 
and characteristic frequencies of CMF are also explicitly obtained. Six CMF characteristic 
frequencies and dispersion strengths against are shown in Fig.2. 

Again, plays an important role in determining the CMF characteristic frequencies 
and the corresponding strengths. Within our model, six CMF characteristic frequencies are 
predicted; three of them are located below lO^Hz and another three do above 10^ Hz, similar 
as that observed in Fig.l. Due to this, the DEP (or ER) spectra are maining characterized 
by two-step rapid changes (or two dominant peaks). For three typical cases, we also find, six 
CMF strengths are all positive in Fig. 2(a); all negative in Fig. 2(b) and three strengths are 
negative and three are positive in Fig.2(c). such behavior will result in co-field or anti-field 
rotation of the particle for ER spectra. We remark that as both the intrinsic dispersion 
and local-field effects are considered simutaneously, we predict four characteristic frequecies, 
characterized by four non-zero dispersion strengths even for = 1/3. 
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III. NUMERICAL RESULTS 



We are now in a position to calculate the dielectric dispersion, and DEP and ER spectra 
based on the model put forward in the previous section. 

In Fig. 3-5, we investigate the effect of depolarization factor on the dielectric response 
[denoted by ee/eo and Oe defined in Eq.(ll)], DEP and ER spectra using Eq.(22) and Eq.(23). 

The dielectric response mainly exhibits two sub-dispersions characterized by two-step 
rapid decrease in Fig.3(a), Fig.4(a) and Fig.5(a). It is evident that the influence of 
(the shape) on the effective dielectric constant (or the effective conductivitty) is significant 
in the low-frequency region where / < IG^Hz (or in the high-frequency region where / > 
IQ^Hz). This is in agreement with previous experimental conclusions that the low- frequency 
subdispersion for e,. depends on the cell shape, whereas the high-frequency is independent 
of it [8]. The dielectric permittivity for oblate spheroidal particles, say = 0.95, in the 
low-frequency region can be as four times as the one for spherical inclusions. Thus, by 
taking into account the shape of cell suspensions, it is possible to reduce some discrepancies 
between the theoretical and experimental results on the dielectric dispersion spectra [8,19]. 

For DEP force, it exhibits strong sensitivity to the depolarization factor (the shape). 
When DEP force is positive (the particles will be attracted towards the field-generating 
electrodes), with increasing L^, such an attractive force becomes weak. However, when the 
DEP force is negative (the particles will be repelled from the electrodes), the larger L^, the 
stronger repulsion will be. Furthermore, we predict one or two shape-dependent crossover 
frequencies at which there is no net force on the cell particle. The crossover frequency is a 
monotonically increasing or decreasing function of L^, dependent on whether the variation 
of DEP force around it is negative or positive [see the inserts in Fig. 3(b), Fig. 4(b) and 
Fig. 5(a)]. This is an interesting result. To the best of our knowledge, the dependence of 
crossover frequency on the spheroidal shape is reported here for the first time. 

For ER spectra, they exhibit two ER characteristic frequencies with both positive 
strengths in Fig. 3(b), both negative strengths in Fig. 4(b), and negative and positive 
strengths in Fig. 5(b) (This has been observed in experimental electrorotational spectra of 
the yeast cell [26,27]). Positive ER torque means the co-field rotation of the particle. In 
this case, with increasing L^, the ER characteristic frequency for ER will be blue-shifted 
accompanied with the decrease of the dispersion strengths. However, when the ER torque is 
negative (the particle exhibits the anti-field rotation) , both the ER characteristic frequencies 
and the strengths decrease concomitantly with L^. Moreover, we find that the variation of 
the geometric parameter L^. does not change the polarity of ER peak. 

We emphasize to point out that, although we analytically predict four and six charac- 
teristic frequencies for dielectric dispersion and ER spectra, these frequencies are mainly 
located in two regions or some of their dispersion strengths are quite small in comparison 
with others, resulting in only two effective characteristic frequencies at which the rate of 
change of the dielectric permittivity (or DEP force) , and the rotational torque peak attain 
their maxima. This can be well understood from Fig.6. 

As we have included the dielectric dispersion form in ei, it would be interesting to 
investigate how the dispersion strength Aei and the intrinsic frequency /i affect the dielectric 
dispersion and ER spectra. 

In Fig. 7, we examine the effect of Aei on the dielectric response (the left panel) and 



7 



ER spectra (the right panel). Clearly, one common fcatTirc of all these dielectric spectra 
is that increasing Aei yields increasing dielectric permittivity in the low-frequency region. 
As previous models do negelect the intrinsic dispersion effect, the theoretical results in the 
low-frequency region are less than experimental reports [8,19]. We think it would be helpful 
to give the closest fit with experimental data by suitable adjustment of both Aei and L^. 
For ER spectra, when Aei = 0, only one dominant peak is predicted as expected. For 
a finite Aei, two dominant peaks at low and high characteristic frequencies arise. With 
increasing Aei, the characteristic low frequency shifts towards a long wavelength; while 
the characteristic high frequency shifts towards a short one, leading to further separation 
between two characteristic frequencies for larger Aei. At the same time, the adjustment of 
Aei can also change the peak value much or less and even dominate the co-field or anti-field 
rotation through changing the polarity of the peak (see the upper right column). 

In Fig. 8, we plot the dielectric dispersion and ER spectra versus the frequency / for 
different /i. It is evident that, /i mainly play a role in the characteristic high frequencies for 
both the effective dielectric dispersion and ER spectra. Generally speaking, the characteristic 
high frequency increases with increasing /i, in accord with the analytical formula from 
Eq.(7). Interestingly, the increase in /i can also lead to the change of the polarity of the ER 
peak in the high-frequency region (see the middle right pannel). 

In order to take into account the mutual interaction bewtween the suspended particles, 
we adopt Eq.(17) instead of Eq.(16). When the mutual interaction is taken into account, 
the ER peak must be reduced, as expected in Fig. 9. With further increasing p, the mutual 
interaction becomes strong, leading to the serious depression of ER peak and the slight blue- 
shift (or red-shift) in the characteristic low frequency (or high frequency). These results are 
also new. 

IV. DISCUSSION AND CONCLUSION 

In this work, we have investigated both analytically and numerically the ac electrokinetic 
behavior, i.e., the dielectric dispersion behavior, the DEP and ER spectra of nonspherical 
cell suspensions with an intrinsic dispersion based on the spectral representation method. 
The dependence of the ac electrokinetic behavior on the depolarization factor is studied in 
detail. We find that an intrinsic dispersion in the dielectric constant of spheroidal cells can 
lead to four and six dispersion in the effective dielectric constant and the CMF respectively. 
We also find that both the intrinsic dispersion strength Aei and the characteristic frequency 
fi can change the polarity of ER torque and thereby causing the co-field or anti-field rota- 
tion. Furthermore, by taking into account the local-field effect from the mutual interaction, 
we examine the ER spectra for various volume fractions, and find that increasing volume 
fractions can result in the decrease of both the strengths and the difference between two 
characteristic low and high frequencies indeed. Thus it is possible to obtain good agreement 
between theoretical predictions and experimental data by the suitable adjustment of both 
the geometric parameters (for example, the particle shape) and the physical parameters (for 
example, the dispersion strength and the characteristic frequency of cell suspensions). On 
the contrary, such a fitting is very useful to obtain the relevant physical information of cell 
suspensions. 

Here a few comments are in order. We obtain four and six characteristic frequencies for 
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the effective dielectric permittivity and CMF. Generally, none of the characteristic frequen- 
cies for the dielectric spectra is equal to those of CMF, as the poles Xi and X2 [Eq.(4)] are 
quite different from the poles Si, S2 and S3 of CMF. In the dilute limit and for spherical in- 
clusions, both Eq.(3) and Eq.(17) yield two characteristic frequencies; the smaller p, the less 
differences between these frequencies [15] are. Numerically, the dielectric dispersion, DEP 
and ER spectra mainly exhibit two rapid changes or two dominant peaks. Here, we also 
show a four-step rapid decrease in dielectric dispersion spectra (sec Fig. 10). However, due 
to small differences between the CMF frequencies or small strengths of some frequencies in 
comparison with others, it is difficult to show all six characteristic frequencies in ER spectra 
clearly. A possible way to achieve this is to consider cell suspensions with large volume 
fractions while being far from the spherical shape. 

We demonstrate theoretically that the shape effect on the effective dielectric constant 
is significant in the low-frequency region, in accord with experimental conclusions [8]. In 
previous work, one of the authors have presented a first-principle study of the dielectric 
dispersion of fission yeast cell suspensions [7] . As the derivation is based on the assumption 
that cell suspensions are in the dilute limit and do not exhibit an intrinsic dispersion, the 
predicted two characteristic frequencies are independent of the volume fractions and some 
discrepancies between the theory and experiment on fission yeast cells still exist. We believe 
that our present formula can be applied also and the discrepancies between theoretical and 
experimental results may further be reduced. 

For shelled spheroidal particles dispersed in the host medium, we have found one peak in 
the electrorotation assay, under the assumption that the ratio of shell to the host dielectric 
constant is real [11]. The theory yields good agreement with the experimental results only 
in the high-frequency region. In fact, for such kind of three component system, the effective 
dielectric response €cs for the coated shells can be firstly obtained in a self-consistent way [21] 
and will exhibit at least one dispersion. The three-component composites can then be 
equivalent to the solid spheroidal particles with ecs including one dispersion term suspended 
in the host medium. Thus our present model can be safely used, and the co-field rotational 
peak in the low- frequency rigion may also be predicted. Finally, we should remark that, for 
nonspherical cells, the intrinsic dispersion may depend on the orientation of a spheroid, which 
will yield dielectric anisotropy in cell suspensions. In this work, we take into account the 
geometric anisotropy as well as the dielectric isotropy. The work including both geometric 
and dielectric anisotropics is in progress and will be reported elsewhere. 
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FIGURES 



FIG. 1. The characteristic frequencies and the dispersion strengths against the depolarization 

factor for p = 0.2 and (a): s = 1.625, t = -0.053, e2 = 78eo, (J2 = lO'^S'/m, Aei = 200eo 
and /i = 1.59 x WHz; (b): s = -0.031, t = 1.25, €2 = 6eo, 02 = 0.025/m, Aei = lOOOeo and 
/i = 2.65 X l{f>Hz; (c): s = 2.78, t = 1.053, €2 = 78eo, ^2 = 2 x IQ-^S/m, Aei = 2500eo and 
/i = 1.9 X lO^ifz. 

FIG. 2. The characteristic frequencies /^i = /ci(si), /62 = /c2(si), /bs = /ci(-S2), /m = fc2{s2), 
fbb = fciiss), fm = /c2(s3),and the dispersion strengths A61 = ViA{si), A&2 = ViS(si), 
A63 = 1^2^(52), A64 = V2B{82), A65 = V^A{s^) and Abe = "^^3^(53) for CMF {b,) versus L^. 
The parameters are same as those in Fig.l. 

FIG. 3. Effective dielectric response (ee/^o and a^) (a) and normalized dielectroplioresis 
(NDEP) and electrorotation (NER) (b) versus the frequency / for the parameters in Fig. 1(a). 
In the insert of Fig.3(b), the crossover frequency f cross is plotted against L^. 

FIG. 4. Same as Fig.3, but for the parameters in Fig. 1(b). 

FIG. 5. Same as Fig.3, but for the parameters in Fig. 1(c). Note that two shape-dependent 
crossover frequencies are predicted. 

FIG. 6. Contributions of dielectric response eei/eo = Aeei/[l+(///ei)^] {i = 1, 4) to the dieletric 
spectra (eg — e//)/eo (a), and contributions of i-th strengths NERi = Abiffu/if^i + *° NER 
(b) versus / for Lz = 0.05. Other parameters are same as those in Fig. 1(c). Note that the effective 
dieletric spectra and the NER spectra mainly exhibit two-step rapid decrease and two dominant 
peaks respectively, although four and six characteristic frequencies are analytically predicted. 

FIG. 7. (ee — e//)/eo and NER versus / for = 0.95 and various Aei. Other parameters are 
same as those in Fig.l. 

FIG. 8. (eg — eij)/eo and NER versus / for = 0.95 and various /i. Other parameters are 
same as those in Fig.l. 

FIG. 9. NER versus / with Eq.(16) and Eq.(17) for = 0.05. Other parameters are same as 
those in Fig.5. Note that Eq.(16) is independent of p. 

FIG. 10. (eg - ei^)/eo and NER versus / for p = 0.2, = 0.95 and Aei = 1 x 10"^ cq. Other 
parameters are same as those in Fig. 1(c). There are four characteristic frequencies, at which the 
dielectric spectra exhibit rapid decrease. 
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